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This article introduces a new multipolar hierarchy for the propagation of the weak-lensing shear, 
convergence, and twist valid in a general spacetime. Our approach is fully covariant and relies on no 
perturbative expansion. We show that the origin of B-modes, in particular on large angular scales, 
is related to deviations of isotropy of the spacetime. Known results assuming a Friedmann-Lemaitre 
background spacetime are naturally recovered. The example of a Bianchi I spacetime illustrates our 
formalism and its implications for future observations are stressed. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

Weak gravitational lensing by the large-scale structure 
of the Universe has now become a major tool of cos- 
mology PQ, used to study questions ranging from the 
distribution of dark matter to tests of general relativ- 
ity [2]. The standard lore [3[ H] states that, in a ho- 
mogeneous and isotropic spacetime, weak lensing effects 
induce a shear field which, to leading order, only con- 
tains £?-modes so that the measured level of 5-modes 
is used as an important sanity check at the end of the 
data processing chain. 5-modes contribution to the ob- 
served shear can be related to intrinsic alignments [5], 
Born correction and lens- lens coupling [6] [7] , and grav- 
itational lensing due to the redshift clustering of source 
galaxies |5]. From an observational point of view, the 
separation of E- and 5-modes requires in principle to 
measure the shear correlation at zero separation [9j [10] 
that can be brought down to the percent-level accuracy, 
e.g. with CFHTLenS data QT]. 

This paper emphasizes that any deviation from local 
spatial isotropy, as assumed in the standard cosmologi- 
cal framework in which the background spacetime is de- 
scribed by a Friedmann-Lemaitre (FL) universe, induces 
5-modes in the shear field. More importantly, and con- 
trary to the above mentioned effects, these 5-modes arise 
on all cosmological scales. Therefore, any bound on their 
level can be used as a constraint on spatial isotropy. This 
is an important signature which, in principle, can be ex- 
ploited in order to disentangle this geometrical origin of 
5-modes from other non-cosmological effects [22]. Since 
it is important for future surveys to predict the level at 
which these cosmological effects produce 5-modes, we in- 
troduce in this work a new multipolar hierarchy for the 
weak-lensing shear, convergence and twist that does not 
assume a specific background geometry. This approach 
will allow us to pinpoint the origin of the 5-modes and, 
in a future work, to access the magnitude of currently 
observed level of 5-modes. 

This work is organized as follows: we start in § [H| 



by reviewing the basic formalism of weak-gravitational 
lensing, which will also help us to set up the basic 
notations and conventions. In 911 Bl we derive the 



evolution equations for the irreducible components of 
the Jacobi map, which are then used to derive the main 
multipole expansion hierarchy in 



IIC 



We then show 
how the standard FL results are recovered (§ III A) 
and dis cuss t he particular case of a Bianchi / (BJ) uni- 
verse (§ IIIB). Finally, we present our conclusion in § IV 



Throughout this paper we work with units in which 
c = h = 1. Spacetime indices are represented by Greek 
letters. Upper case Latin indices such as {/, J, K, ... } 
vary from 1 to 3 and represent spatial coordinates. Fur- 
thermore components of vectors on a spatial triad (a 
set of three orthogonal spatial vectors which are nor- 
malized to unity) are denoted with lower case Latin in- 
dices {z,j, fc, ...}, whereas the screen projected (two- 
dimensional) components are represented by indices 
{a, 6, c, . . . } which vary from 1 to 2. 



II. MULTIPOLAR HIERARCHY FOR 
WEAK-LENSING 

A. Description of the geodesic bundle 

A crucial quantity for weak-lensing is the electromag- 
netic wave-vector, k^ = d^w, where w is the phase of the 
wave. In the eikonal approximation, k^ is a null vector 
(k^kjj, = 0) satisfying a geodesic equation (k v V v k^ = 0). 
Moreover, if we assume that V M V ' v w = V v V^w for any 
scalar function w (torsion- free hypothesis), it follows that 
its integral curves x^(v) defined by k^(v) = dx^/dv, 
where v is the affine parameter along a given geodesic, 
are irrotational (V^fc^j = 0). Second, we consider a fam- 
ily of null (light-like) geodesies collectively characterized 
by x^(v,s), where s labels each member of the family. 
We adopt the convention according to which v = at 
the observer and increases toward the source. There is a 



wave- vector for each geodesic, that is k^(v : s) = dx^/dv, 
and the separation between the geodesies is encompassed 
by the vector r]^ = dx^/ds connecting two neighbor 
geodesies (see Fig. (fil)). Hence, we first derive the dy- 
namics for a reference geodesic, and then the dynamics 
for the deviation vector. 




FIG. 1: Representation of two null geodesies of the light bun- 
dle, iff is the projection of iff in the plane spanned by the 
basis {e a }. The dotted curve represents the worldline of the 
observer comoving with u^ . The geodesic bundle is thin so 
that its transverse dimension has not been depicted and it 
converges at the observer. 



We suppose that the light-rays converge to a funda- 
mental observer comoving with the four-velocity u^ of 



matter, which is normalized such that u 
observer measures a redshift z given by 



z(v) 



_ (k^) v 



so that the energy of the incoming photon is 



-1. This 



(1) 



(2) 



In this work we adopt the perspective of a photon going 
to the past, which means that in a local Lorentz frame, 
where u^ = ( — 1,0,0,0), we have k° = dt/dv = —U. In- 
cidentally, this suggests that we introduce of a "reduced 
wave-vector" through 



fc" 



dx» 
dv 



U' 1 ^ 



(3) 



in order to simplify our expressions x 



1 Since d0 = —dt, the new parameter v is simply the negative of 
the proper time t, reflecting our choice of perspective in which 
the observer sheds light on the source. 



At each position x M of a given geodesic we can asso- 
ciate a direction vector n whose components are n M , and 
defined from the reduced wave- vector through 



with 



fc" 



u^n^ = , 



n^n* 1 



1. 



(4) 



(5) 



At the observer, n° = n(v = 0) is the spacelike vector 
pointing along the line of sight 2 . However, since we now 
have 



dv 

dv 



U, 



it follows that we can either choose (n°,v) or (n°,v) as 
independent set of variables to parameterize the geodesic, 
which correspond to two slices of the past lightcone. As 
we shall see below, the use of v simplifies the derivation 
of the multipolar expansion for the weak-lensing observ- 
ables. 

At a given point of the geodesic, it is necessary to add 
two vectors to u and n in order to obtain a complete basis 
of the tangent space. We choose these two vectors n a , 
with a = {1,2}, to be orthonormalized and orthogonal 
to u and n, that is they are defined by 



<™b/* 



5ab 



<*V = <1 



= 0, (a = 1,2). (6) 



Since n and n a comprise a three-dimensional orthonor- 
mal basis, we can simplify the notation by defining 713 = 
n so that we can collectively write n^ = {n^} i= i^^. Note 
that at the observer we can again define n° = rii(v = 0) 
with a remaining rotation freedom around n° for the 
choice of n° a . 

We now introduce the screen projector tensor 



^flU 9fIU ~T~ ^[L^V TlflTlw; 



(7) 



which projects any tensor on the two-dimensional surface 
orthogonal to the line of sight. Thanks to the orthogo- 
nality relations (|6| , the basis can be parallel transported 
along the null geodesic as [12] 



S^WVvnl = 0. 



(8) 



At v = 0, each n° of the geodesic bundle can be as- 
sociated to a spherical basis and this can be used to 
fix the rotational freedom. Indeed, for each n° there 
will be a unique choice of n°(n°) and n^n ) if we set 
{e£,eg,e°} = {n° ',n?,n§}. The integration of Eq. ^ 
then allows to define this basis at each point on the past 
lightcone, i.e. to determine 71^(71°, £), or, equivalently, 



2 Note that our definition of n^ differ by a minus sign from that 
ofRef. [13] 



{e r , e<9, e^}(n°, v) everywhere. This prescription empha- 
sizes the importance of introducing a reference triad as a 
way of identifying these projection effects; see Fig. [2j 

At this point it is convenient to introduce the helicity 
basis defined as 



1 



1 



e± = n± = — = (e e T ie^) = -y= (rii =F in 2 ) • (9) 



Their components in the n a basis read simply 
n a ± = n ± .n a = ±={5$T^) 



(10) 



and are, by construction, constant. 

We now note that any event on the lightcone is 
uniquely specified by (n ,^), i.e. it is of the form 
^(n ,^). This means that any local quantity X(x^) 
evaluated on the lightcone can be seen as a function 
X(n°,v). The redshift defined in Eq. ([I]) is also a func- 
tion of (n ,^), and U propagates as (see e.g. Ref. [13 ) 



dlnC/ 

dv 



H«(n°,v) 



(11) 



where the parallel Hubble expansion rate along the line 
of sight is defined by 



is the screen projected Riemann tensor. The linearity of 
Eq. ( 15 ) implies that 



^>) = w(trL' 



where the Jacobi map U a b satisfies the Sachs equation [3, 

m 



subject to the initial conditions 

dV a 
W)=0, ^(0) = ^. 



(17) 



(18) 



In order to proceed, we need to decompose both V^ 
and 1Z a b in their irreducible pieces. We start by decom- 
posing the projected Ricci tensor into a trace and a trace- 
less part as 



K ab = U 2 (Mat + W ab ) 



(19) 



where 1Z and W a b are related to the Ricci (R^ u ) and Weyl 
(Cupav) tensors through: 



H«(n°,v) = WVnU,,. 



(12) 



Using the standard 1 + 3 decomposition of V ^u vi it takes 
the general form 

H\\(n°,v) = \& + (Ji lv n»n v + A^, (13) 

where 6, a^ and A M are the expansion, shear and accel- 
eration of the flow vJ 1 . All these quantities are evaluated 
on [x^(n°,v)] and are thus functions of n(n°,v) on the 
past lightcone. 



Ti=--R^k\ W ab = C^ au kek°ny^ (20) 



and where 



■Lab — ^fiiy n a n b 



(21) 



is the identity matrix of the screen space. Note again 
that W a 6, as well as 1Z and 7Z a b, are evaluated on the 
central geodesic and thus W a b [^ /i (n°,'0)] = yV a b(n°,v). 
In terms of the electric and magnetic parts of the Weyl 
tensor, given respectively by [20] 



B. Shear, twist and convergence propagation 

The purpose of this section is to derive an equation 
governing the shear, twist and convergence of a light- 
ray bundle without specifying the spacetime structure. 
The evolution of the deviation vector n^ is given by the 
geodesic deviation equation 



dV 
~d^~ 



R* A „aek»k a rf , 



(14) 



where R^ U a(3 is the Riemann tensor. This equation can 
be rewritten in terms of its component on the screen basis 
{n a } as [3] 



where 



d<2 Va -, b 

^2 = KabV , 



i^ab — RfivapK nl ^^5 



(15) 



(16) 



the projected tensor W a b becomes 

Wab(n°, v) = -2n£<> [£^ + B^e av (n)] , x . 



(n° , v) 



(23) 
In the expression above, () stands for the traceless part 
with respect to I ab . e^^n) is the antisymmetric tensor 
in the projected space and is defined as 



e/xi/(n) = vPefo 



(24) 



Now, W a b being a spin-2 field, it can be decomposed in 
the helicity basis d9| as 



W ab (n°,v) = -2^2 W X (n°,v)n^. 



(25) 



A=± 



This decomposition emphasizes once more that the two 
components W x are functions of (n°,v) alone, because 



\G r , e#, &tp\o 



<~&ch 



^fejj°^d 



*c e 




tr ?£d~~-^ 







^\m:^> 



Background 



|e r , GQ^ Qip\(jlO^ 



FIG. 2: Any position on the past lightcone can be considered as x M (n°,£). While quantities such as S^ u are local, quantities 
such as Wab depend on the local basis at x^(n°,v) via the projection on n^n^y, see Eq. (23). Observational quantities are 
however defined in terms of n° so that one needs to relate the basis {e r , e#, e v } in (n°, v) ana in v = 0. The relation n(n° , v) 
induces "projection effects" and a non-local relation between quantities like Sim and W a 6- Once a background spacetime 
is chosen, its symmetries simplify the comparison. For instance, a B/ spacetime provides a natural triad of Killing vectors 
associated to its principal axis. One can use this "global reference" to relate the local S 2 in x^(n°,v) to the observer's S 2 by 
comparing them in the reference S 2 . 



they are evaluated on the lightcone. Recall that the n^ 
are constant so that we can use either n or n° in Eq. (25 ). 



We now decompose the Jacobi map in terms of a con- 
vergence ft, a twist V and a traceless shear j^ as 



V n 



Klab + Ve a 



' lab i 



(26) 



where 



tab 



2iTl [a n o] 



All these quantities are defined on our past lightcone so 
that we can also think of them as functions of (n°,v). 
The shear, being also a spin-2 field, is naturally decom- 
posed similarly as 



jab(n°,v)= ^7 A K,^X' 



(27) 



A=± 



Finally, by inserting the decompositions ( 25p7 ) in the 
Sachs equation (17) we find the desired equation of evo- 
lution 

nx- 7 +) \ 

Wt- 7 +] -(28) 
W ± (k±iV)J 

Note that, in practice, the integration of this system re- 
quires the evaluation of the past lightcone structure in or- 
der to determine rii(n ,v) and then H\\{n°,v), 7Z(n°,v) 
and W ± (n°,v). 




C. Multipole expansion 



Equation (28) is composed of scalars (ft, V, 1Z and H\\) 
and spin-2 fields (7 ± and W ± ) defined on the sphere. 
The former can be naturally decomposed in a basis of 
spherical harmonics as 



K(n°,v) = ^2 K£m(v)Y im (n°) 

£,m 

V(n°,v) = J2 v em(t>)Yem(n°) 
e,m 

■R(n°,v) = Y, n tm{b)Ye m (n°) 
e,m 

H n (n°,v) = J2 h ^(v)yem(n°) 



(29) 
(30) 
(31) 
(32) 



i.m 



The latter, being spin-2 fields on the sphere, can be ex- 
panded on a basis of spin-weighted spherical harmon- 
ics [15] as 

W ± (n°,v) = ^2[S £m (v)±iB £m (v)}Y^(n°),(33) 



£,m 



w±/ 



£,m 



Note that E- modes are those having parity ( — 1)^ while 
5-modes have parity (— 1)^ +1 [16]. 

It is important to keep in mind that we are adopting 
an observer-based point of view so that all quantities are 
expressed in terms of (n°,v). In general, n(n°,v) ^ n°, 
with the obvious exception of e.g. FL spacetimes or for 
an observer at the center of symmetry of a Lemaitre- 
Tolman spacetime. Part of the difficulty is thus contained 
in the determination of these coefficients, which include 



projection effects from the geodesic structure. 



When inserting these decompositions in Eq. ( 28 ) , prod- 



ucts of spherical harmonics will appear on the r.h.s. They 
can be simplified using standard relations between spin- 
weighted spherical harmonics (see Appendix [A|. It fol- 
lows that, in terms of multipoles, the equations of evolu- 
tion for the convergence, twist and shear take the follow- 
ing general form 



d Eirr 

dv 2 



dv 2 



d 2 Kerr 

dv 2 



Qrnrnirri2 ljZe imi ~ hl imi — I [5 L E£ 2 m 2 -\- lS L B£ 2 m 2 ) — 2K,£ irni (5 L £t 2 m 2 +^L^2 

+2Ve imi (-iSlSe 2m2 + 5 L Be 2 m 2 )] 

%imi — feimiTT j (S L B£ 2 m 2 ~ iS L E£ 2 m 2 ) ~ 2K,£ irni (5 L B(, 2 m 2 ~ l^L £ 



■{' 



2 y^mm\m 2 

Ly u 1 e 2 
2Ve imi (5 L iBe 2m2 + 5 L £i 2 m 2 )] 



livn 2 ) 



s^tmm\m 2 



i^l\m\^l 2 m 2 i^iirni 



dne 2 






2(— l) mi C i2 ™ 2 £i mmi [S L (Ei irni Ei 2 m 2 + Bi irni Bt 2 m 2 ) + i^L {E^ irni 8i 2 m 2 — Ei irni Bi 2 m 2 )\ } 



(35) 



(36) 



(37) 



d 2 V e „ 
dv 2 



J s~imm\m 2 



'R-£ 1 m 1 V£ 2 m 2 ~ feimi 



dV (2 



dv 



+ 2( — l) mi Cg™g 2 £l Tnrni [5 L i(E£ 1 m 1 ££ 2 m 2 + Bl 1 m 1 Bl 2 m 2 ) ~ S L {Bt irni & 2 m 2 ~ E£ 1 m 1 B£ 2 m 2 )\} 



(38) 



r 



where 



<^ = 



[l±(-l) L ]/2, L = £^£^£ 2 (39) 



and an implied sum over £i, £ 2 , mi, and m<2 is under- 
stood. This multipolar hierarchy for weak lensing, which 
does not rely on a particular background spacetime - 
and on any perturbative expansion - has never been de- 
rived before and sets the basis for general studies of the 
constraints on anisotropy and inhomogeneity from the 
weak- lensing 5-modes. 

As soon as the spacetime has a non-vanishing Weyl 
tensor, E- and 5-modes are generated due to the cou- 
pling of the Weyl tensor to the convergence and twist. It 
shares some similarities with the Boltzmann hierarchy for 
the cosmic microwave background (see e.g. Refs. [HJ[T7]) 
but one needs to keep in mind that 7v!^ m , /i£ m , £^ m , Bi m 
are non-local quantities since they have to be evaluated 
on the geodesic. 



III. APPLICATIONS TO SPATIALLY 
HOMOGENEOUS UNIVERSES 

A. Standard FL case 

In order to illustrate the formalism we consider the 
standard case of (flat) FL spacetime with linear pertur- 



bations. At the background level, the metric of the FL 
spacetime takes the simple form 



ds 2 = -dt 2 + a 2 (t)S IJ dx I dx J 



(40) 



This spacetime enjoys 3 translational Killing vectors 
i e i}ie{x,y,z} = {-^r}ie{x,y,z} which define everywhere 
a natural Cartesian basis. By normalizing these vectors, 
we can then define a triad of vectors e^ whose compo- 
nents are e/ = Sf /a (and their associated 1-forms e % 
whose components are e l I = <5}a), that is a set of three 
orthonormal space-like vectors (and forms) that can be 
used as a global Euclidian basis. The set of vectors n?, 
which was a priori only defined at the observer's position 
can then be defined everywhere by imposing that their 



components n, 



03 



n? .e J in this reference basis remain 



the same everywhere. This enables to compare rii(n ', v) 
to n° even though these sets of vectors are defined first 
at two different points of spacetime, as illustrated in the 
right part of Fig. |2| 

At the background level, the Weyl tensor vanishes (i.e. 
£\iv — an d S/xiv = are at least of order 1 in perturba- 
tions) and the Ricci scalar, VS°\ depends only on time. 
For this spacetime n(n°,v) = n° for all v so that the 
only nonzero multipolar coefficient h^ is the monopole 



do) 

^00 



H 



(41) 



where the dot refers to derivative with respect to t. From 
the expression above and the fact that dv = — dt, we find 
from Eq. (11 ) that U oc a -1 . It then follows from Eq. ^ 
the well known result 1 + z = ao/a. Moreover, since 
Bi°2 = 0, it follows from Eqs. (I37I and S that 



Indeed, first order vector and tensor modes would gener- 
ate 5-modes since then Br^ ^ 0. 



The equation (37) for the convergence has r.h.s. 



(o) 



S -0 

Kq and Vq satisfy the same second order homogeneous 

equation of the form 



n t 



(0) ,(0) 



00 



00 



dv 



.(i) 



^ (1 V 0) 



(51) 



as usual, since the other terms in Eq. (37) are at least 



d 2 X 



(o) 
00 



dv 2 



HY (0) 
/^oo^oo - H d ~ 



,(0) 



,(0) 



where XL stands for either 



v oo 



JO) 
^00 



or V r 



(0) 



00 



(42) 



The initial 

.(0) 



of second order in the perturbations. For the twist, the 
argument is similar but TJLlVqq = and the initial con- 



.(i) 



ditions (18) imply V}^ = 0. So, in conclusion 



conditions (18) then lead to a homogeneous ft 00 , given 
by the usual angular distance, and a vanishing twist so 
that 






7^0 



v (l) - 



0. 



(52) 



JO) 

""00 



D A , 



V r 



(0) 



00 



0. 



Then, one concludes that 



E, 



(0) 



Bi°2 = 0. 



(43) 



(44) 



At higher order, £ a & and B a b are non- vanishing (note 
that one cannot simply drop out Bt m in the hierarchy, 
even for pure scalar modes), which leads to I?-modes as 
well as twist. Moreover, projection effects and couplings 
induced by hg m need to be included; see Ref. [18] for the 
case of second-order perturbations. 

The absence of 5-modes at first order in perturbations 
are due to the fact that 



At first order in the perturbations, the perturbed met- 
ric with only scalar perturbation reads in the Newton 
gauge 



1. B 



(i) 



for scalar modes, 



2. at this order we can work in the Born approxima- 
tion. 



ds 2 = -(1 + 2&)dt 2 + a 2 (t)(l — 2^)5udx I dx J , (45) This latter point is extremely important since otherwise 



where $ and \I/ are the two Bardeen potentials. The 
projected Ricci tensor is of the form [12] 



even if B^ = the dependence n(n°, v) would generate a 



it 



(i) 

ab 



-D a D b ($ + *) 



(46) 



where D a is the covariant derivative on the 2-sphere. It 
follows that 



B 



ab 



0. 



(47) 



In the Born approximation (i.e. n(n°, v) = n°), only 
Hqq ^ so that the r.h.s of Eqs. (35JJ36) involves only 



2 / om0m2 



non- vanishing B a b 0- Indeed, in Eqs. ( 35p7 ) part of the 
difficulty lies in the determination of the coefficients 7£g m , 
him, £im and Be m that depend on the whole geodesic 
structure, as we shall now illustrate. 



B. Example of a Bianchi / 

We now consider the case of a spatially homogeneous 
but anisotropic universe described by a Bianchi / space- 
time for which the metric takes the form 






£0£ 2 



. Thus £2 — f- an d L is even (i.e. S L =0). As 



a conclusion, in Eq. ( 35 ) for the propagation of the E- 



ds 2 



-dt 2 + a 2 (t)-f IJ (t)dx I dx-' 



(53) 



modes, the only remaining term on the r.h.s. is 



n: 



(0) _^(0)_ 



00 



6 00 



dv 



?W 



2k (0 V (1) 

Z ^00 c £m 



(48) 



where the coordinates have been chosen so as to diago- 
nalize Ju(t). This solution is spatially homogeneous and 
the spatial shear 



while in Eq. (36) for the propagation of the 5-modes it 



cru 



ld7/j 
2 dt 



(54) 



n t 



(0) 
00 



<(0) 
TO 



dv 



B 



(i) 

£m- 



(49) 



So we see that only i^-modes are sourced, while 5-modes 
would need to be initially non-zero to be non-vanishing 
today, 



*£!*o 



B, 



(i) 



0. 



(50) 



characterizes the spatial anisotropy, a 2 (t)jjj being the 
spatial metric, a(t) is the volume averaged scale factor 
and 9 = 3H = 3d /a (see Refs. [19 for notations and 



3 The "standard" convergence and shear are —k^/k^ and 
'ab l K ^ * n our notations. 



properties). It follows that the kinematical quantities 



entering Hn in Eq. (13) are 



6 = 3tf, a 7J ^0, 



0. 



(55) 



Similarly to the FL case, 
Killing vectors {e/}/ G{a . >3/ ^} = 



this spacetime enjoys 3 



idx 1 * 



le{x,y,z} 



which de- 



fine everywhere a natural Cartesian basis. Normalizing 
these vectors, we can then also define a triad of vectors 
Si that can be used as a global Euclidian basis. And sim- 
ilarly to what has been done in the FL case, the set of 
vectors n° can then be defined everywhere by imposing 
that their components in this reference basis e$ remain 
the same, in order to allow the comparison of rii(n ,v) 
ton?. 

However, contrary to the FL case, one has to consider 
(i) the non- vanishing background electric Weyl tensor, 



£f) = Hau + -a 2 lu - a IK a^ (56) 

while the magnetic part is identically null, 



B 



(0) 

u 



0, 



(57) 



and (ii) the fact that at background level rii ^ n° (un- 
less in the particular case of geodesies along one of the 
three proper axis), which induces projection effects so 
that/i^O. 

The triad ni{n° , v) is related to the reference triad n° 
by a rotation defined by three Euler angles as 

ni(n°,v) = fl/K/^K, 

= Jn g (7)/Jn§(/3)/Jn g (aO fe << (58) 

where the Euler angles are also functions of (n°, v). The 
determination of a, /3 and 7 requires the integration of 
the geodesic equation in the Bianchi spacetime. 

Then, for a typical tensor T^ v at an event x M , 
its projection orthogonally to n, i.e. its components 
T ± [x /i ,ni] = T ± [x /i (n°,'0)] in the helicity basis n±, can 
be related to its projection at the same event x M , but 
orthogonally to n° with components T^[x^,n^] in the 
helicity basis n° ± . For a spin s tensor, this transforma- 
tion reads (see details in appendix [B| in general 

T^x^m] = exp(±i^)exp(/3 a D a )T ± (^,<), (59) 



with 



a + 7 and (3 = /3[nJ cos 7 + n° 2 sin 7] . 



For a homogeneous spacetime, the dependence in x^ of 
T^- reduces to a time dependence. Eq. (59) evaluated for 
a rank-2 tensor (that is s = 2) is needed to account for 
the projection effects in the definition of W ± [x /i (n°, v)]. 
Similarly, Eq. ( 59 ) in the case s = (that is for a scalar 



field) is needed for the projection effects of H\\ [x^(n° ', v)\ 



and 1Z[x fJ, (n°,v)]. The Weyl tensor having only a non- 
vanishing electric part (with only non-vanishing compo- 



nents £ x 
one has 



£ yy and £ zz in the natural Cartesian basis), 



W?(v,n°) 



ra=0,±2 



£°2mHW^{<) (60) 



with 



-2 



CO 

c 2,±2 




^yy) 



C11 



(61) 
(62) 



The projection of the electric Weyl tensor has a direc- 
tional dependence for £ = 2 and m = 0, ±2. However, 
the directional dependence of <p and f3 a in Eq. (59), 
i.e. the projection effects, sources and mixes E and 
B modes at higher £, as for CMB polarization E/B 
modes mixing [21 . This projection effect also induces 
non- vanishing 1Zi m terms even if the background Ricci 
is homogeneous. 

To go further and understand how this mixing of E 
and B modes arises, let us assume that <r 2 /0 2 is small, 
so that we can work at first order on this parameter (we 
can think of B/ has a homogeneous perturbation of FL). 
Then, the geodesic equation and the parallel transport of 
n a [Eq. d8k] lead to 



dv 



n 



S ik a kj n j 



&il 



dv 



nt 







(63) 



and thus at lowest order one easily obtains that <\> ~ 
and /3 a (n°,v) ~ J* D a a(n°,v')dv' . Here a(n,v) = 
oij(v\nJn J j2 can be thought as a lensing potential and 
Eq. d59b for W gives 



W ± [x^n z ] ~ [l + /3°D a ]W^(r7,n?), 



(64) 



similar to the form for linearized lensing in FL [17j [25] on 
light polarization. a(n°, v) obviously contains only £ = 2 
multipoles. Because of the derivative coupling, using [T7] 



D V ]~) a \ r ^ zS \ ^ T ±s s~immim2 \r±s 

U a Y ixrnx U Y i2m2 - 2_j L Mil* ^££^2 Y £m 



(65) 



where 



Lu^ = 5 Mh + 1) +4(4 + 1) - £(£ + 1)] m 



and further defining 



±s jmm\m2 



s jmm\m2 r ±s sv 

J- pp. 0- = L1U1Z2 ^ 



±s /^iram\m2 



one can convince oneself that, at background level, terms 
such as 



°£m — C tm "•" 1 U 1 £ 2 



£ ae imi dv'\ dtS^ m2 (67) 



and 



Appendix A: Spin-weighted spherical harmonics 



B 



(0) 



2 jmm\m2 



-l Z l 



<r £imi dv' S L £i 2m2 (68) 



are expected when extracting the E and B modes out of 
Eq. 64 . Thus a multipolar £ = 4 B-mode will appear. 



One needs however to rely on the full transformation ( 59 ) 



so that the E- and B-modes shall be generated for larger 
€s. Similar sources arise from 1Z and H\\ for which pro- 
jection effects will generate non- vanishing TZ^ and h^. 
A full analysis, including perturbations and magnitude 
estimations will be presented in Ref . [22] . Our argument 
sketches the expected effects that arise from the higher 
multipoles induced by the background Weyl tensor and 
the fact that n ^ n°, an effect that cannot be neglected 
even in the Born approximation for anisotropic spaces. 
Besides, in BJ spacetimes, the amplitude of vectors and 
tensors is of order of the shear times the amplitude of the 
scalars, another source of 5-modes. 



IV. CONCLUSION 

We have provided a new multipolar hierarchy for weak 
lensing. Our formalism, which is fully covariant, does not 
rely on perturbation theory nor on the choice of a back- 
ground spacetime. It allows us to relate the property 
of the shear to symmetry properties of the background 
spacetime and discuss the generation of £> -modes. We 
have argued that a violation of local isotropy is expected 
to leave a B-mode signature on all scales. This result 
is important for future surveys, such as the Euclid mis- 
sion [23] (early results on the B-modes have already been 
obtained from CFHTLS [26] and DLS [24 and we can 
forecast that Euclid will typically decrease the error bars 
on the 5-modes by a factor of order 10-40 on scales rang- 
ing up to 40 degree, that is in the linear regime where 
astrophysical sources of 5-modes are expected to be neg- 
ligible) and may us allow to set new constraints on the 
deviation from spatial isotropy on cosmological scales. 
The quantitative computation of the level of 5-modes 
expected on large scales, where the gravitational dynam- 
ics can be considered linear, for a Bianchi universe is 
currently being investigated [22] and requires to study in 
details the cosmological perturbation theory beyond the 
analysis of a scalar field [19] . 
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We gather here a few important identities and relations 
between spin-weighted spherical harmonics used in this 
text. The reader is referred to Ref. [15] for more details. 

Spin- weighted spherical harmonics form a complete set 
of orthonormal functions on the sphere, satisfying 



Jd 2 nY* ±s (n) em Y^,(n) 



SfPt $n 



(Al) 



An important identity, used in particular to derive Eqs. 
(|35l-(|38|), is 






±s 

£m 



(A2) 



where 



±S/^mmim2 



S /T 



±S 



d*nY*r(n)Y timi (n)Y*'(n). (A3) 



Using the transformation of spherical harmonics under 
parity, it can be shown that the above coefficients satisfy 



=ps s^mm\m2 / i\L ±s s-imm\m2 

Ly U 1 £ 2 ~ v L ) ^£ii£ 2 

where L = £ + £ x + £ 2 . 



(A4) 



Appendix B: Expansion of tensors on the sphere 



The transformation ( 58 ) of the triad is easily computed 



for the tangent space basis when the helicity vectors are 
used. Indeed we first note that for a general rotation 
around an axis n, the corresponding helicity basis at the 
point of the sphere, that is n±{n) 1 transforms as 



Jn(<p) n± 



exp Tl( ^ n± , 



(Bl) 



i.e., it is a diagonal matrix in this basis. Furthermore, 
when the axis of the rotation does not coincide with the 
direction defining the helicity basis, the rotation is still 
diagonal in the sense that, if we consider a general rota- 
tion J not necessarily around the axis n, then 



nt(J • n) • J • n±(n) = 



(B2) 



or J • n±(n) ex n±(J • n), where we recall that n±(J • 
n) is the helicity basis at the point of the sphere which 
is the image of n by J. It comes essentially from the 
fact that rotations defined as 50(3) are not only keeping 
orthogonality conditions, but also orientations of triads 
as their determinant is required to be 1. 

In order to use these interesting properties, we refor- 
mulate the transformation (58), which is written with 



rotations around the axis of a fixed frame, by its form 
where the rotations are performed around the axis of the 
rotating frame. In that case it reads 

m(n°, v) = U(a)/ J n 'AP)i k J n -M£ n° (B3) 



with n 2 ° = J n ° (7) • ^2 an d n 3° — Jn'°^) ' n s- Using the 



transformation rule (Bl) and the property (B2), we de- 
duce that the transformation rule for the helicity vectors 
is just 



n±(n°,v) = exp^ a+7 ) J n > (P) ■ n°± . 



(B4) 



It then proves convenient to express a rotation as a 
parallel transport. Indeed under any rotation of angle if 
around an axis n ro t, a tensorial quantity T in the tan- 
gent space at a point n equat0 r of the corresponding equa- 
tor (that is such that n eqU ator * ^rot = 0) is transformed 
exactly as if it were parallel transported with the vector 
^transport = ^™rot x n eqU ator. Let us note this parallel 
transport 7^ tra ort (T) . We insist that this rephrasing 
of a rotation as a parallel transport is valid only on the 
equator of the rotation. In our cas e, this is enough to 
reformulate our transformation (B4) as 



with 



n±(n°,v) = exp^ a +^ Tp(n° ± ) 



(3{n°) = /3[n?cos7 + n£sin7] 



(B5) 



(B6) 



and 7/3 being the parallel t rans port along (3. 

Let us apply the result (B5) to obtain an expression 



for the expansion of a tensor on the sphere. For a rank- s 
tensor T, its projection orthogonally to n defines a tensor 
field on the sphere. For instance, the rank-2 tensor W^ 



defines a tensor field on the sphere S^S^W pa given that 
the screen projector S pv depends on the position n on 
the sphere of directions. In the evaluation of the geodesic 
deviation equation, we are led to express the components 
T ± = T ± [n(n°, v)] of a tensor field at a point n(n°, v) in 
the helicity basis n±(n° , v) in function of its components 
T^- = T^(n°) at a reference point n° in the helicity basis 
n±. This expansion is obtained as follows 



T • Tlzp . . . fl^\n{ri ,v) 

Tp 1 (T • n T . . . n T )\ n o 



(B7) 



n^ . . . n^ 



= exp^*^^ [exp ((3 a D a )T\ n o] . n T ...n° T 

From the first to the second line, we have used that a 
scalar field (the components of T) evaluated in n or its 
parallel transport back along (3 evaluated in n are equal. 
From the secon d to the third line, we have used the trans- 
formation rule (B5) of the triad. From the third to the 



fourth line, we have used the exponentiation of the paral- 
lel transport in terms of covariant derivatives D a on the 
2-sphere. Then, with a common abuse of notation (see 
for instance the discussion at the end of Ref. [25]), which 
we also consistently use in Eq. (65), this is rewritten in 
a short form as 



T ± = exp 



±is(o;+7) 



exp(/rD a )T d 



(B8) 
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